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The Magnetospheric MultiScale Mission (MMS) is scheduled to launch in late 2014. Its 
primary goal is to discover the fundamental plasma physics processes of reconnection in the 
Earth’s magnetosphere. Each of the four MMS spacecraft is spin-stabilized at a nominal 
rate of 3 RPM. Traditional spin-stabilized spacecraft have used a number of separate modes 
to control nutation, spin rate, and precession. To reduce the number of modes and simplify 
operations, the Delta-H control mode is designed to accomplish nutation control, spin rate 
control, and precession control simultaneously. A nonlinear design technique, Lyapunov’s 
method, is used to design the Delta-H control mode. A global spin rate controller selected as 
the baseline controller for MMS, proved to be insufficient due to an ambiguity in the 
attitude. Lyapunov’s design method was used to solve this ambiguity, resulting in a 
controller that meets the design goals. Simulation results show the advantage of the pointing 
and rate controller for maneuvers larger than 90 deg and provide insight into the 
performance of this controller. 


Nomenclature 


B 

nj 

= 

angular velocity of the body, B, with respect to the inertial frame, N 

COq 

= 

desired spin rate magnitude, constant 

[■/] 

= 

3x3 Inertia matrix 

B 

X 

= 

applied external torque to spacecraft in the body frame, B 


= 

angular momentum of spacecraft in body frame, B 


= 

desired spacecraft angular momentum in body frame, B 



angular momentum error in the body frame, B 

Ho 


desired spacecraft angular momentum magnitude 

[I] 

= 

3x3 identity matrix 

X 

= 

maximum eigenvalue of [J] 

s 

= 

desired major axis pointing direction in inertial space, N, fixed in inertial space 

s 


instantaneous unit vector representing the desired major axis pointing direction in inertial 
space resolved in the body frame, B 

^ 1:3 

= 

scalar positive constant 


= 

scalar positive constant 


= 

desired body pointing vector, fixed in body frame, B 


= 

controller error vector in the body frame, B 


= 

external applied torque in the body frame, B 


I. Introduction 

T he Magnetospheric MultiScale Mission (MMS), being designed and built by NASA’s Goddard Space Flight 
Center (GSFC) in Greenbelt, MD, is scheduled to launch in late 2014. The primary goal of MMS is to 
discover the fundamental plasma physics processes of reconnection in the Earth’s magnetosphere. Each of 
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the four MMS spacecraft will be spin-stabilized at a nominal spin rate of 3 RPM and have the same complement of 
in-situ science instruments. Science instruments are placed at the end of two 15 m long axial coilable booms, two 5 
m long radial booms, and four 60 m long axial wire booms. All spacecraft will be launched together on the Atlas V 
launch vehicle. The attitude control system (ACS) for this spin-stabilized spacecraft is not as simple as traditional 
spin-stabilized spacecraft. MMS's Delta-H control mode is the primary attitude mode and is responsible for 
simultaneous nutation control, spin rate control, and precession control; it’s design is the main topic of this paper. 
Section n gives an overview of the attitude control system (ACS). Section III delves into the baseline Delta-H 
control design, presents a correction to the baseline controller, and discusses thruster-based implementation. Section 
rV shows simulation results comparing the baseline and augmented controllers and give the reader a better feel for 
the Delta-H mode performance. A conclusion is presented at the end of the paper. 

II. Attitude Control System 


A. ACS Design 

The design philosophy is to have the ACS as simple as possible. There are two main control modes, one for 
attitude maneuvers (Delta-H mode) and one for orbit adjustment maneuvers (Delta- V). There is a Sun Acquisition 
mode, which is only used in the event of an off-nominal separation. A complete open-loop mode is offered for 
thruster checkout and spin-up during boom deployments. No control actions are performed in Science mode, the 
default start-up mode. Figure 1 is the ACS control mode transition diagram showing all the ACS control modes. 



Figure 1. MMS ACS Control Mode Diagram 


B. ACS Components 

There are three ACS sensors on each MMS spacecraft. A digital sun sensor (DSS) is used as the primary Sun 
Acquisition attitude sensor, as well as for attitude and rate verification during Checkout mode operations. There is a 
4-headed star sensor that provides up to 4 attitude solutions to the Delta-H and Delta- V control modes. An 
accelerometer measures translation during Delta- V maneuvers. There is no angular velocity sensor; an angular 
velocity is derived from a Kalman Filter. 

MMS has twelve thrusters used for attitude and orbit control. Eight 4 lb thrusters are directed radially, while four 
1 lb thrusters are directed axially. Refer to Fig. 2 for a thruster placement diagram. 

- - INSERT Figure 2. Thruster placement on MMS — 

C. Attitude Maneuver Operations Concept 

After sun acquisition, the ACS is required to maintain the spacecraft’s spin rate at 3 RPM while maneuvering up 
to 180 degrees to attain mission attitude. At mission attitude, the ACS must adjust the spacecraft spin rate from 3 to 
7 RPM to accommodate the boom deployment sequence. At the conclusion of the deployment sequence, the 
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spacecraft spin rate is back to 3 RPM. Then, during science operations, the ACS must precess the spacecraft up to 
two degrees every two weeks, maintain a constant spin rate, and minimize the core body nutation. Traditional spin- 
stabilized spacecraft have used a number of separate control modes to control nutation, spin rate and precession. To 
reduce the number of modes and simplify operations, we seek a controller which can perform all three 
simultaneously. The remainder of this paper details the MMS Delta-H controller design. 

ni. Delta-H Controller Design 

Spacecraft attitude dynamics are governed by Euler's Equations: 

[JJ + 02 X[J\ 02 = X , (1) 

where is the spacecraft angular velocity, [J] is the spacecraft inertia matrix, and is the sum of 

external torques on the spacecraft. Theoretically, if is small then Eq. 1 can be linearized. Looking at the 

attitude maneuver operations concept for MMS, there is a requirement to: maneuver 2 deg, spin up the spacecraft 
from 3 to 7 RPM, and maneuver up to 180 deg. Equation 1 can be linearized for the first and second maneuver. For 
the first maneuver, linearizing about a constant spin rate, and neglecting higher order terms, leads to a decoupling of 
the spin and transverse dynamics. For the second maneuver, assuming small transverse rates leads to a single 
integrator plant for the spin dynamics. The third maneuver cannot be linearized. Since we are seeking one controller 
to perform all attitude maneuvers, a nonlinear controller will be designed. Many options exist for nonlinear 
controllers. A Lyapunov design approach is chosen. 

Conceptually, we want the Lyapunov function to be proportional to an energy-like term. If we ensure the change 
in the Lyapunov function is consistently decreasing, we can claim the controller drives the system to a minimum 
energy state. There are varying levels of nonlinear stability and the most desired is Global Exponential Stability 
(GES). GES means that from any initial state, the system will reach the single stable equilibrium point exponentially 
fast. See Ref 1 for more details on nonlinear design strategies. 

A. Global Spin Rate Controller 

Reference 2 presents a globally stable Lyapunov controller for spin stabilized spacecraft. The Lyapunov based 
controller is designed for simultaneous spin, nutation and precession control. Let a candidate Lyapunov function be 
a continuously differentiable function defined in a domain Z)c |R^ that contains the origin (as presented in Ref 
2 ): 


{’H-H/sf 


2 — 

X|J| 

-1 

[/]] 


\ 





( 2 ) 


where, 

^K=[j] is the spacecraft angular momentum vector, is the spacecraft angular velocity, ^H.j) 

is the desired angular momentum, is the error between the current and desired angular momentum, [ / ] is 

the identity matrix, [j] is a 3x3 inertia matrix, A is the maximum eigenvalue of [j] , Hq is a constant 
representing the desired angular momentum magnitude, [K^] is a symmetric positive definite gain matrix, and 
is an instantaneous unit vector representing the desired major axis pointing direction in inertial space resolved 
in the body frame. Because desired major axis pointing direction in inertial space resolved in the inertial 

frame, N, is fixed in inertial space, ^ J varies by: 


B^ N B^^B^ 

S-— (JO X s 


(3) 


If the derivative of the Lyapunov function is negative then the Lyapunov function will tend towards zero, i.e., a 
minimum energy state. Therefore, to prove Lyapunov stability, we must show that V approaches zero. 

Taking the derivative of Eq. 2 yields: 
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V{x) = ^Hl^HyH^[K,YH , (4) 

Substituting Euler's equation yields: 

Simplification leads to: 

V=("RyR\K^'^^\ ( 6 ) 

Further substitution and simplification follows: 


= Yr-^Rd+^R\Kp'^^^t ( 7 ) 

Substituting for : 

=('££-'i£„+'i£|x[^r'-|/|]r'i (8) 

Expanding: 

=(-'iii,+'fl:xkl''r"i (9) 

Substituting for and leads to: 


=(-///j+[j]vxkrT"i ■ (>•» 

Recall is a unit vector representing the major axis; therefore Hq—\{Dq . Substitution and simplification 
yields: 


= — cOq^S) ( 11 ) 

F=X(Vf'l 02) 

If the external torque applied to the spacecraft is — ,i e., in the opposite direction of the rate 

error, and [ K^\ is a 3x3 symmetric positive definite gain matrix, then: 

F=-x(vr[xjv 03) 

Equation 13 is negative semi-defmite; it must be negative definite to ensure GES by Lyapunov's Theorem. 
Therefore, we seek bounds on Eq. 13 such that we can guarantee it to be negative definite. From the definition of the 
induced 2 -norm and the fact that the eigenvalues of a symmetric matrix are real, we get the following bound: 

F<-XA.„.([/:,])||V||:, x,.([/:j )>0 (H) 

Defining the term 0^3 = A A,^^^ ([ ^c] ) substituting leads to: 
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, N B 2 

V <-«3 O), 2' 


V«3>0, V"^o)felR' 


(15) 


Equation 15 is negative definite for all values of > which implies this controller is GES by Lyapunov's 

Theorem. 

If the torque applied is aligned with the rate error, V is negative, and the final state will be V = 0. At V=0, the 
following occurs: 


— — \ oJq^s (16) 

Equation 16 shows that the desired momentum is achieved. In addition: 

^i£^[A[jr'-[/]]^i£=A'tt)o'v[x[jr'-[/]i^5 . (17) 

Equation 17 shows that a spin about a major axis is achieved. 

However, an ambiguity exists in the direction of major axis. If the major axis is the z-axis, and you want to point 
it at > then the controller can end up pointing the body at either +z or -z. This manifests itself in the controller 
development. Recall the rate error is given by = 02 — COq 5* , which indeed is driven to zero. However, 

inspection of the rate error and recalling Eq. 3 reveals that ^ J is kinematically coupled with the spacecraft’s 
attitude, i.e., ^ ^ s , where is the direction cosine matrix from the inertial 

frame to the body frame. Since the problem fonnulation for this controller did not include an attitude term in the 
Lyapunov function, it cannot guarantee convergence of the attitude. The controller still remains GES in rate-space 
(i.e., the rate error is driven to 0 from any initial condition), however it is not global with respect to attitude-space. 
Because MMS is concerned with attitude (pointing) as well as rate, an additional term must be added to Eq. 2 to 
account for the ambiguity. 

B. Global Pointing and Spin Rate Controller 

Augment the Lyapunov function, Eq. 2, with an attitude-dependent term, while maintaining the same conditions 
as Eq. 2: 


Tenn 1 


Tennl 




(18) 


Tenn 3 


where , and k 2 are scalar gains whose sum must equal 1, and is the desired body pointing vector 

(fixed in body-space). 

Term land Term 3 in Eq. 18 are components of the controller in Section III. A. Term 2 is the additional term 
added to resolve the attitude ambiguity. By adding Term 2, which is proportional to energy with respect to the 
attitude, we have selected a Lyapanov function that may guarantee rate control, nutation control, and precession 
control as well as guarantee the rate and major axis are pointing in the desired direction. 

Taking the derivative of Eq. 18 and simplifying (similar to the procedure taken in Section III. A.) yields: 


V = \(ki + k2] 


N B P 

O) — COq 


\r 


k^-\rk2 


T — k^ e T 


target vector 


( 19 ) 
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where the target vector is a weighted sum of the body pointing direction and desired spin direction, 
k^ — X[k^-\-k^ , and is the difference between the current spacecraft angular velocity and the weighted 
target vector. If we let — , where is a 3x3 symmetric positive definite gain matrix, then: 

( 20 ) 

From the definition of the induced 2-norm and the fact that the eigenvalues of a symmetric matrix are real, we can 
bound Eq. 20: 


Defining the term — k^ ( [^c] ) substituting leads to: 

V<-a,\Y£„ V«3>0, V^gGlR' (22) 

Equation 22 is negative definite for all values of , which implies this controller is GES by Lyapunov's 
Theorem. If the torque applied is aligned with the error vector, V is negative, and the final state will be V = 0. 
From Eq. 18, when V = 0, the desired momentum, spin about a major axis, and the desired pointing direction is 
achieved. 

Looking at the error vector, , reveals that the controller is essentially a proportional+derivative (PD) 
controller, since the error vector is composed of a term proportional to attitude and a term proportional to the rate. 
Even though a nonlinear design technique is used, the end result is a PD controller, a common result of nonlinear 
control design. 


rV. Simulation Results 

A. Advantage of Global Pointing and Rate Controller 

Simulation results reveal the attitude ambiguity in the rate controller. When we perform a maneuver greater than 
90 degrees, the spacecrafts is shown to despin and spin up in the opposite direction. Comparing the rate controller 
with the pointing and rate controller, shows that adding the attitude dependency resolves the ambiguity in the rate 
controller. 

- INSERT Figure x. Spacecraft rate for two controllers. - 

B. Performance 

The Delta-H controller performs well. Figure xx shows a 60 deg precession maneuver on the unit body sphere. 
Notice that the angular momentum vector goes along a straight line from the initial to the desired state. Also notice 
that the body z-axis spirals towards the final state. The spiraling action is indicative of induced nutation from the 
execution of the precession maneuver. Once the final state is reached, the angular momentum and the body z-axis 
align and come to rest, meaning that the controller was able to control both precession and nutation. 
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Figure xx.60 deg maneuver showing angular momentum on body 

sphere 

Equation 1 8 gives two scalar gains, and k 2 . Notice that if =1 and =0, then the controller 

collapses to the global rate controller, and all that it entails. The ratio of these gains represents the weighting 
between pointing vs. nutation. The target vector is a linear combination of the major axis pointing direction, , 
and the body pointing direction, . Weighting more heavily on results in better pointing accuracy and 

weighting more on results in better nutation control. Therefore, these gains should be selected appropriately 
based on the desired goal of a class of maneuvers. 

V. Conclusions 

Traditional spin-stabilized spacecraft have used a number of separate control modes to control nutation, spin rate 
and precession. To reduce the number of modes, and simplify operations, a controller to perform all three was 
designed. The Delta-H control mode accomplishes nutation control, spin rate control, and precession control 
simultaneously. In designing the controller following Lyapunov's technique, we found that a global Lyapunov rate 
controller was insufficient for our application because it did not guarantee global stability of the attitude. Therefore, 
a pointing and rate controller was developed using Lyapunov re-design that proved to guarantee the desired 
objective. Simulation results show the benefit of the pointing and rate controller for maneuvers greater than 90 deg. 
Simulation also shows the general performance characteristics of this controller and provides some insight into 
selecting the controller gains. 

Although theoretically complex, this controller is essentially a PD controller. This GES controller should be used 
as a baseline for future spinner missions, as the design is theoretically sounds, and its implementation is quite 
simple. Future work involves significant flight software testing in preparation for flight in 2014. 
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